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We study, in some detail, the linear stability of closed timelike curves in the Godel universe. We 
show that these curves are stable. We present a simple extension (deformation) of the Godel metric 
that contains a class of closed timelike curves similar to the ones associated to the original metric. 
This extension correspond to the addition of matter whose energy-momentum tensor is analyzed. 
We find the conditions to have matter that satisfies the usual energy conditions. We study the 
stability of closed timelike curves in the presence of usual matter as well as in the presence of exotic 
matter (matter that does satisfy the above mentioned conditions). We find that the closed timelike 
curves in the Godel universe with or without the inclusion of regular or exotic matter are stable 
under linear perturbations. We also find a sort of structural stability. 



o 
o 

(N 



oo 

(N 



X 



PACS numbers : 04.20.Gz, 04.20.Dw, 040.20 Jb 



> . I. INTRODUCTION 

o 

The Godel universe, a spacetime in which the matter takes the form of a pressure-free perfect fluid with a negative 
, cosmological constant, is the most celebrated solution of Einstein field equations that contains closed timelike curves 
(CTCs). This spacetime has a five-dimensional group of isometries which is transitive. The matter everywhere rotates 
| relative to the compass of inertia with the angular velocity proportional to the square root of the matter density fl[ 0] ■ 
Dynamical conditions for time traveling in this spacetime are not sufficient to exclude the existence of CTCs |3(. 
*q | The Godel metric has some qualitative features like the projections of geodesies onto the 2-surface (r, d>) being simple 
q-i closed curves. This property can be extended to a set of metrics called Godel-type. It is possible to show explicitly 
that when the characteristic vector field that defines a Godel-type metric is also a Killing vector we have closed 
timelike or null curves [J] . The study of the geodesic motion of free test particles in these Godel-type spacetimes can 
be extended to a family of homogeneous Godel-type spacetimes Q. Giirses et al. [6| introduce and used Godel- types 
metrics to find charged dust solutions to the Einstein field equations in D dimensions. In the Godel spacetime timelike 
geodesies are not closed as was independently proved by Kunt Q and Chandrashekhar and Wrigth 8] . However there 
^ ' exist timelike curves subject to an external force that are closed. 

Exact solutions of Einstein-Maxwell equations that contain CTCs at least for some values of the parameters are 
studied in Q. It turns out that magnetic fields can give rise to non-trivial chronology violations. A sufficiently 
large magnetic field can always ensure chronology violation. The spacetime of an infinite rotating cylindrical shell of 
charged perfect fluid contains CTCs [l(|. For the conditions for the existence of CTCs in the spacetime associated 
to a rigidly rotating cylinder of charged dust see Ivanov fill . In spacetimes with conic singularities that represent 
cosmic strings we can have CTCs for one spinning string jl2| | and for two parallel moving strings (l3| . 

The stability of Godel's cosmological model with respect to scalar, vector, and tensor perturbation modes using a 
gauge covariant formalism is considered in [l4j]. It is found that the background vortical energy contributes to the 
gravitational pull of matter, while gradients add to the pressure support. The balance between these two agents 
effectively determines the stability of Godel universe against matter aggregations. 

The possibility to have no violation of causality for geodesies that in finite intervals of time goes back in time 
(dt/ds < 0) was considered by [l5| . As it was already mentioned, the existence of CTCs contradicts the usual notion 
of causality. Beyond the usual paradoxes, it seems to induce physical impossibilities, like the necessity to work with 
negative energy densities. One could speculate that these impossibilities will be eliminated by quantum-gravitational 
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effects. All our experience seems to indicate that the physical laws do not allow the appearance of CTCs. This is 
that, essentially, says the Chronology Protection Conjecture (CPC) proposed by Hawking in 1992 [la ]. 

Other possible explanation for CTCs is that the metrics studied until now are unrealistic [13] ■ But, there exists an 
exact, asymptotically flat solution of vacuum Einstein equations which contains CTCs that can represent the exterior 
of a spinning rod of finite length [181 ] . 

The existence of CTCs, in principle, should be a matter of experimentation. If the General Relativity predicts them 
in a physically reasonable situation and they are not found, we will have that this theory is in trouble. If they are 
found we will have a bigger problem, our usual notion of causality will need a deep revision. 

The possibility that a spacetime associated to a realistic model of matter may contain CTCs leads us to ask how 
permanent is the existence of these curves. Perhaps, one may rule out the CTCs by simple considerations about their 
linear stability. Otherwise, if these curves are stable under linear perturbations the conceptual problem associated to 
their existence is enhanced. 

In this first paper about linear stability of CTCs we consider the Godel case due to its historical relevance (first 
paradigm for CTCs) as well as its mathematical simplicity. We intend to study, in subsequents works, the stability 
of CTCs in Godel- type metrics 0,0], Bonnor metric [HI and other relevant spacetimes. All these metrics presents 
their own peculiarities that deserve a special treatment. 

The paper is divided as follows, in Section 2, we review in some detail the CTCs in the Godel metric. In particular 
we present the CTCs and its corresponding forces in the usual Godel coordinates as well as in Cartesian coordinates. 
In Section 3 we study the stability of CTCs when they suffer a small perturbation. We solve the linear system of 
equations for the evolution of the perturbation in the original Godel coordinates in an exact form. For future reference 
the same system in Cartesian coordinates is solved numerically for a representative range of initial conditions. In 
Section 4 we present a simple extension (deformation) of the Godel metric that corresponds to the addition of matter. 
We look for the conditions for geodesies to be timelike and to be closed. We prove that these two conditions can not 
be satisfied at the same time. We present the CTCs in an explicit form together with their associated forces. We find 
that the addition of matter either usual or exotic does not change the stability of the Godel CTCs. In the last section 
we summarize and discuss some of our results. 



II. GODEL METRIC, GEODESICS AND CTCS 

In this section we review the CTCs in Godel metric following mainly references [l| 0] 0]. In standard coordinates, 
X 1 * = [t, r, <p, z], the Godel metric is [H, 



dP - dr 2 + ( sinh 4 r - sinh 2 r) dip 2 + 1\FL sinh 2 rdtdp -dz 2 , (1) 



where (3 is a constant that describes the vorticity of the four-vector u M = [1,0,0,0], — oo < i < oo, — oo < z < oo, 
r > and <p £ [— it, ir}. The limits ip — ~tt and <p = ir are topologically identified. 
The metric ([TJ satisfies the Einstein field equations, 

R^v — + A-<7jui/ = —Tfj, v , (2) 

with T^ v = pu^u", A = - ^ and p = /3 2 . We use units such that c = 8nG = 1. 
Let us denote by 7 the particular closed curve given in its parametric form by, 

i=i , r = r , <pe[-7r,7r], 2 = 0, (3) 

where to and ro are constants. When 7 is parametrized with an arbitrary variable u the condition to be timelike is 
^JiT^r > 0- The curve 7 is timelike when ro > log(v / 2 + 1). The four-acceleration of this CTC is 

= 5% sinhr coshr (2 cosh 2 r - 3) <p 2 . (4) 
In the "Cartesian" coordinates X^ — (i, x, y, z) defined by 

y = (cosh(2r) + sinh (2r) cos p) , (5) 
fix = (sinh(2r) simp) y, (6) 

<p + ( ^r_2f )=arctan / e -2r tan <M (?) 

2 v 2V2 V 2 / ' w 

z = z, (8) 
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with — oo < t < oo, — oo < x < oo, — oo < z < oo, and oo > j/ > we find that the line element (JTJ) reduces to 



ds 2 



dt 



V2dx 



-i 2 



dx 2 + dy 2 

(Py) 2 



dz 2 



In these coordinates the geodesic equation gives us, 



•• 2 . . v/2 . . 
t = -ty+ xy, 

y py 
x = -0V2 iy, 

y = (3V2ix+ X , 

y 

z = 0, 



(9) 



(10) 

(11) 
(12) 
(13) 



where the overdot indicates derivation with respect to the proper time s. This system of equations admits the first 
integrals, 



t 



2y-yo 

V2C ' 



0y 1 \ 
x = -givo - y), 

y = -g{x-x ), 
d 

z= c> 

where xq, yo, d and C are integration constants. The solution of this system is 



V2 



2 arctan (VN tan a — 



N+l 



2VN 



+ t 



2 V TV 7] sin a cos a 

9 , AT • 2 ^ X °1 

cos a + N sin cr 

cos 2 tr + TV sin 2 f7 ' 
2cr d 



/3 y/y'o - V 2 



with 



N = 



2C 

yo-y 
yo + y 



- so), 



where to is an integration constant, and r\ is given by r/ 2 := (x — xq) 2 + (y — yo) 2 - 
The relation X^X^ = 1 (with z = 0) reduces to 



x 2 + y 2 



1. 



From (fIl |) -([T6 ]l we get 



Therefore 



y 2 = 2(C 2 +7 ? 2 ). 



(14) 

(15) 
(16) 
(17) 

(18) 

(19) 
(20) 
(21) 

(22) 
(23) 



(24) 

(25) 
(26) 
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The condition for a geodesic to be closed, i.e., £(— 5 ) = i(§ ), x(— 5 ) = £(5-) e y(— 5 ) = y(f ), is 



= N + 1 =*> y'o = r,v 2 

that contradicts (|26[) . Hence we have no timelike closed geodesies in the Godel universe. 

In order to obtain a CTC in Cartesian coordinates we use another differential equation for t obtained from 
and equations (O and (|20]l. Then 



(27) 



whose solution is 



t = 



2V2 



■ V2, 

t + — (yo - v) 



c 2 + v 2 
c 2 



arctan (^N tan 



y Q V2 - VC 2 +?7 2 



+ t . 



These timelike curves are closed when 



C 2 = 2(y Q -^/y 2 - V 2 ) 2 - V 2 



(28) 



(29) 



(30) 



By replacing these parametric equations for the CTC into the geodesic equations we find that the four-acceleration 
satisfies the relations 



a° = -(x-x )— — , a 1 = Xy(x - x ), a 2 = -Xy(y - y), a 3 = 0, 
P 



(31) 



where 



X - ^2 



a/2(C2+7, 2 )-j/ 



(32) 



The right hand side of (|3"Tj) can be interpreted as the components of a specific external force associated to 7. 

It is instructive to compare the force given in these two systems of coordinates. In Godel standard coordinates we 
have g vv (f> 2 = 1. Therefore along the CTC 7 we have 



U 2 



4( sinh 4 ro — sinh 2 ro) 



From equation ([4]) we can write the non zero component of the force in following way, 

- 1 (3 2 sinh ro cosh ro (2 cosh 2 ro — 3) (3 2 sinh 2ro ( cosh 2ro — 2) 
F = - 



4(sinh 4 r - sinh 2 r ) 2((cosh2r - 2) 2 - 1 

The non zero components of force in Cartesian coordinates are, 

2V2 



F" 



VP 



(x-x )F\ 



F 1 = ^-(x-x )F\ 



F 

F 3 



2jL (yo -y)F\ 



= 0. 



On 7 the constants in ([22l and (|23|) are now N = e 2r °, yo = cosh2ro and 77 = sinh2ro. Therefore 



F 1 



P 2 v(yo - 2 ) 

2((2/ - 2) 2 - 1) 



(33) 

(34) 

(35) 
(36) 

(37) 
(38) 

(39) 
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By using the condition (130)) we can cast A as 



Hence 



that is force calculated in (T3Tj) . 



X - (2/0-2)2-1- (40) 



F 1 = | A, (41) 



III. LINEAR PERTURBATION OF CTCS IN THE GODEL UNIVERSE 

A generic CTC 7 satisfies the system of equations given by 

X fi =F»(X), (42) 

as 

where D J-J is the covariant derivative of the vector field () Q along j(s) and is a given external force. 

Let 7 be the curve obtained from 7 after a perturbation £, i.e.jA^ = X* 1 Let e Q be a given basis. In this basis 

(|42| is represented by the equation S± = F, where ^ = (^ + r^u /3 -u^)e a , F = F a e a , u = u Q e Q and u a = X a . In 
order to find the behavior of £ we calculate the variation, in first approximation, of both sides of = F. We find, 

— = + 2r£ i WV + 1^ Vu")e a + F a Se a (43) 

<5F = i^e a + F a £e a . (44) 

where = Jpc- Comparing the last two equations we get the system of differential equation satisfied by the 
perturbation £ 

^ + 2r^«WV + r^ iA £V^ = F^. (45) 

This last equation can be cast in a manifestly covariant form by noticing that its left-hand-side is the well known 
geodesic deviation equation as pointed out in [l|| and its right-hand-side is the Lie derivative of the force along £ M . 
In standard coordinates X^ = [i, r, <p, z], the system (|4"5|) reduces to 

^ +ai 1 = 0, 



where 



(/3 sinh 2 r a/2) / (cosh ro y sinh 2 r — 1 



b = /3V2COS11 r /y sinh 2 r — 1, 

c = /3 cosh r (2 cosh 2 r — 3)/\J sinh 2 r — 1 , 

d = (3 2 cosh 2 r (2 cosh 2 r - 3) 2 /(2 sinh 2 r (sinh 2 r - l) 2 ), 

e = /3/(sinhr cosh r y^sinh 2 r — 1) (47) 

The solution of system of equations is, 

£° = -a(c 3 sin(ws + 04)/^ + ts) + ci s + c 5 , 

= c 3 cos(ws + c 4 ) + r, , 4g . 
i 2 — — e(c 3 sin(ws + c 4 )/w + rs) + c 2 s + c 6 , 
£ 3 = c 7 s + c 8 , 
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where Ci, i = 1, ...,8 are integration constants, oj — yjd — ab — ce, and r = —(feci + cc 2 )/lu 2 . In order that the 
perturbed curve, 7, remains on the plane z = we take initial conditions such that c-j = c§ = 0, i.e., £ 3 = 0. The 
solution (|48|) shows the typical behavior for stability, we have vibrational modes untangled with translational ones 
that can be eliminated by a suitable choice of the initial conditions. 

For future reference, we study the linear stability of the CTCs in Cartesian coordinates. In these coordinates the 
system (|45[) can be written as 



where (ay) and (bij) are given by 



Kj) = Q 







1- a 02 £ 2 + 




- b Q2 e 




+ an^ 1 ~ 


^a 12 i 2 + 




- b 12 e 


e = a 20 i° 




1- a 22 £, 2 + 




- b 22 e 


= 0, 










2(3Se~ 2ra 
D 




V2S 




[1 + (y - 2)De 2r ° 


y/2(3 2 Se- 
D 2 


4r 










2L> 



,-2r 



D 



(49) 



/%o - 1) 



2£Se_ 

IT 



-2r 



(50) 



(&«) 



_A\/2 
e" 2ro A 



— iy/2S 



m 

c 2 

2e- 2r « S A 



cos 2 a + e~ 4r ° sin a 

A 



-4r. 



2e- 2r o 



'0 sin (7 



sm ct 



2/0 



,2 fl2 



c 2 



(51) 



and 



D = cos 2 a 



-4r 



sin a, S — r/ sin er cos cr. 



As before, we take £ 3 = in order that 7 does not leave the plane z = 0. 

To describe the behavior of perturbed curve 7 we introduce two distance and one angle functions, 

R\ = (x~x ) 2 + (y-y ) 2 , 

R 2 = {t-t ) 2 + {x-x ) 2 + (y-y ) 2 , 

y-vo 



arctan ■ 



x — Xq 



(52) 



(53) 
(54) 

(55) 



The first function is a constant equal to 77 when x and y are on the CTC 7. The second represent a "radius" in 
spacetime and the third one is an angle on the usual space. The variation of these functions along 7(5) are, 



sr 2 = [(x-x )e + (y-vo)e}/v, 

sr 3 = [(t-t )f + (x-x )e + (y~yo)e]/R3, 

S<p = [(x - xo)^ - (y - ycO^A?. 



(56) 
(57) 
(58) 

To study these functions we solve the system (|49[) by running the independent variable cr, first from —ir/2 to and 
second from ir/2 to 0, we recall that the points —ir/2 and ir/2 on the curve are identified. We keep the initial position 
of the perturbation equal to zero and analyze the behavior of SR 2 ,S(j) and SR3. The perturbation initial velocity is 
taken each time with only one component different from zero. 

In Fig. [T]we show graphics that represent the variation of radius R 2 of perturbed curve, with ro = 1 .5(77 ~ 10) and 
(3 = 1. For the first graph the initial conditions are: a) £ M (-§ ) = [0,0,0,0], ^(-f) = [10^ 3 , 0,0,0] (solid line), b) 



FIG. 1: Variation of the radius R2 of the perturbed curve, with ro = 1.5 and (3 = 1. The initial conditions for the first graph are: 
a) ^(-|) = [0,0,0,0], ^(-f) = [10~ 3 , 0,0,0] (solid line), b) = [0,0,0,0], = [0, 10" 3 ,0,0] (dot-dashed line), and 

c ) £"(-§) = [0,0,0,0], £ M (-§) = [0,0, 10~ 3 ,0] (dashed line). And initial conditions for the second graph are: a) £ M (§) = [0,0,0,0], 
= [-10-3,0,0,0] (solid line), b) the ^(f) = [0,0,0,0], ^(f) = [0, -10~ 3 , 0, 0] (dot-dashed line), and p(f) = [0,0,0,0], 
^M(|) = [0,0,-10- 3 ,0] (dashed line). 



^(-f) = [0,0,0,0], C M (-f) = [0,10" 3 ,0,0] (dot-dashed line), and c) £^(-§) = [0, 0, 0, 0], ^(-f ) = [0,0,10" 3 ,0] 
(dashed line). And for the second are: a) ^(f) = [0,0,0,0], ^(§) = [-10~ 3 , 0, 0, 0] (solid line), b) the = 
[0,0,0,0], ^(f) = [0,-10~ 3 ,0,0] (dot-dashed line), and = [0, 0, 0, 0], <^(§ ) = [0,0,-10^,0] (dashed line). 

In Figs. [2]and[3]we show graphics that represent S(f> and SR3, respecticely, computed with the same initial conditions 
used for 5B,2. 



x10~" 6 * x10~" 6 * 




-4 1 ' 1 1 1 1 1 1 1 -2 1 1 1 1 1 1 1 1 1 

-1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



FIG. 2: Variation of the angle <j> for the perturbed curve with the same initial conditions of Fig. [T] 



We consider a generic perturbation 8u small when > < 0.01, i.e., (5u) 2 « 8u. All the variations presented in 
the graphics satisfy by large this smallness conditions. We have typical vibrational modes for 8R 2 and SR3 and mainly 
translational modes for 5<f>. We tested this smallness condition for a significant variety of initial conditions and values 
of ro, we find that this condition is always satisfied even for curves that are near the bound given by tq > log(%/2 + 1). 



IV. CTCS AND MATTER 



The material source of the Godel metric is a pressure-free perfect fluid. In this section we study the persistence 
and stability of the CTCs when we add matter with pressure or tension. 
Let us consider the metric, 

ds 2 = (dt + ^dxf - j^(dx 2 + dy 2 ) - dz 2 (59) 

that for h\ = hi = 1 reduces to Godel metric. This metric can be considered as a deformation of the original Godel 
metric (See Appendix A). 
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FIG. 3: Variation of R3 with same initial conditions of Fig. [T] 



First, to understand the physical meaning of the changes introduced in (|59p we compute the associated energy- 
momentum tensor from the Einstein field equations with the same cosmological constant as in the Godel metric 
(A = -/3 2 /2). We get, 



I3 2 (Zh\-2hl+h%) V2f3hi(2hl-hl) 
t3 2 {hi-h\) 



2h* 





2/4 









f3 2 (h$-hl) 









2 {hl-2hl+hj) 
2h| 



(60) 



By solving the eigenvector equation, 



(61) 



we find the eigenvalues: Aq 



A(3 2 
B 



Ai — At — 



C£ 
B 



2 D3 2 
and A 3 = where A = 3/if - 2/i| + hj,B = 2h\ , C = h%-h\ 

B 



and D = h 2 — 2h\ + h\. The timelike eigenvector u M = [1,0,0,0] is associated to Aq, and the spacelike eigenvectors 



X» 



hi h 2 



(3y 

= [0,0, —,0], and Z^ = [0,0,0,1] are associated to Ai, A2, and A3, respectively. We 

"2 



can write the energy-tensor in its canonical form as, 

= XqU^u" + X^X^X" + Y»Y V ) + X 3 (Z»Z"). 



(62) 



In order to have realistic matter, the eigenvalue Ao that represents energy density, and the eigenvalues Ai = A2 
and A3 that describe pressures (or tensions) are restricted. Ao must be non-negative (weak energy condition) and 

Ao — Ai — A2 — A3 > (strong energy condition) 0]. The weak energy condition is satisfied when h\ > 2,l2 3 ,l2 . If, 

furthermore, we have h 2 > -f the strong energy condition is also satisfied. 

Simple modifications as the one presented in (|59p of the Godel metric written in standard coordinates are not 
associated with an energy-momentum tensor with a simple physical interpretations nor to geodesic equations that 
can exactly be solved. This justify our use of Cartesian coordinates. 

Returning to the analysis of the curves in this new geometry, we have that the geodesic equations are, 



2h\ ■ . s/2h x 



(2hl - h\)xy = 



PV2hi • . 2h{- h\ 



h 2 2 



-tx 



yh 2 



( X ) 2 --(y) 2 =0 



(63) 

(64) 

(65) 
(66) 
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By appropriated choice of the constants of integration, the first integrals to these geodesic equations can be written 
as ([Tt |) -([T7 1) replacing i by 

i=-^-[2yhl-y (2hl-hl)]. (67) 



Therefore 



t = -V 2arctan(ViVtan a)h x - - — \ -- ]=J-a +^o- (68) 

/? /ii 2-//V 



We have timelike geodesies for 



that are closed whenever 



yo = cosh 2,o < 7 == r (69) 



2/i? 
2h\ - hi 



y = cosh 2r = 1 (70) 



Hence, as before, we have no closed timelike geodesies in the spacetime whose metric is given by (|59[) . In this metric 
we can find equations for CTCs in the same way as before. For z = 0, we have 



2^2 



The condition for closed curves is 



The components of force are, 



arctan I v TV tan a ) h\ — — 7=^=^= — = — 

C 2 = 2hl(y Q -^yl- 1 ff-hW. 

F^ = -( X - XQ ) ^f^ , 

Fp = (x-x ) y \{h 1 ,h 2 ), 
F p = -V Hh 1 ,h 2 )(y Q -y), 
F* = 0, 



+ No- 



where 



X(hi,h 2 ) 



C 2 h\ 



We have CTCs when C 2 > 0. From (l?2"j) we find the relation between h\ and h 2 to have CTCs, 



hi > 



V 2 hl 



2(yo - 1) 2 



(71) 
(72) 

(73) 

(74) 
(75) 
(76) 

(77) 
(78) 



From the equations for the CTCs it is possible to analyze its stability under linear perturbation. As before, we 
write the coordinates of the perturbed curve as X^ — X^ + and find a system like ljl5|). where (dy), (fry) are now 
given by 



(aij) 



Ch 2 



2V2f3 2 Se- ir Qh 1 
D 



2V2 shim 



4/3Se- 2r °N 2 

w 



V20 2 e- 2r Oh 1 (y o D-e- 2r Q) 2/3( e - 2r ° (h 1 -N 2 )-D(h 2 1 ~N 1 y ) 
W D 



2V2h 1 (2h 1 -N 1 + (y N 1 -2h 1 )De 2r °) 
2 

-2/3(e- 2 '-Q(jV 2 +/ H )-g(fe 1 +Af 2 y )) 
D 



2l3Se- 



D 



(79) 
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(bij) = 



Xy/2 



e- 2ro (iVi - hi) + D{ht - N iy 0) 



e- 2r " A 
D 





Dh 2 C 2 



~ 2r ° S 
D 



48 2 N9 e~ 2ro 



2e~ 2r ° 



cos 2 a + e 4, sin 2 a 



h\C 2 v D 
P 2 



h 2 C 2 



3 -2ro 



Ni(—-yQy + 



4h 2 S 2 e~ 4r °' 
D 2 



(80) 



Ni = 2h\ - h 2 , N 2 = h\ - h\, and A is like in ([77|. As before, | 3 = and we do £ 3 = 0. 



i i i 
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FIG. 4: Similar curves as in the Fig[T]with the same initial conditions. But, now we have ordinary matter with parameters h\ = 
and /12 = 




FIG. 5: Similar curves as in the Fig. [T]with the same initial conditions. But, now we have exotic matter with paramcrtcr h± = V2.7 and 
h 2 = V2A. 

The analysis of the linear stability of CTCs in this new universe is made in the same way as in the Godel's case. 
Several initial conditions are tested and the results for usual as well as for exotic matter are very similar to those 
obtained before. We tested the stability of CTCs in various scenarios. In Fig. [4] we plot the function SR2 for 
the same initial conditions used in Fig. [1] and hi — vTT, hi = \J\.\. We see a similar behavior as in the pure 
Godel case. In Fig. Owe show 8R2 for the same initial conditions for exotic matter (Ao — \\ — A2 — A3 < 0) with 
h\ = \/2.7, = \/2.4. We see that the linear stability is kept. 

For some values of hi and /12 the condition C 2 > [see, Eq. ([72")) ] is not true. Therefore it is possible to have 
spacetimes with a metric like (|59j) with no CTCs for ro > log(l + V2) and any kind of matter (ordinary or exotic). 

In Fig. [S] three lines divide the (ft. 2 ,, ft 2 )-plane in four regions of interest that are marked with /, //, 177 and IV. In 
the region above the dot-dashed line [h 2 = (2h\ — /i|)/3] we have matter with positive energy density. Above the solid 
line [h\ = {h^)/2] the strong energy condition is fulfilled. Under the solid line we have exotic matter. The dashed 
line (straight line) separates the values of hi and /12 for spacetimes with CTCs (above) to the ones that do not have 
CTCs (bellow). This line is h\ = 0.862030830483155/i| that is obtained replacing r = 1 in (75]). 
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Therefore, when {h\,h\) £ I the spacetime contains CTCs and ordinary matter. For (h\,h\) £ II the spacetime 
contains ordinary matter, but not CTCs. If hf) € III the spacetime contains CTCs and exotic matter. And when 
(h^hi) £ IV the spacetime contains exotic matter, but not CTCs. The isolated point represent the Godel universe 
(hijhz) = (1, 1). As we can see in this figure we have an open neighborhood of the point (1, 1) (Godel spacetime) 
where the matter is ordinary, the CTCs are present and are stable as depicted in Fig. [6j 



3.5 




-0.5 1 ' 1 1 1 

0.5 1 „ 1.5 2 2.5 

FIG. 6: For values of (h^, h\) S / we have spacetimes with CTCs and ordinary matter. When (h^hf) £ II the spacetime contains 
ordinary matter, but not CTCs. If 6 /// the spacetime contains CTCs and exotic matter. And for £ IV the spacetime 

contains exotic matter, but not CTCs. The isolated point represent the Godel universe (hi, I12) = (1, 1). 

Loosely speaking, we have structural stability of a vector field when the equations that it satisfies are slightly 
changed we also have a small change in the trajectories represented by the vector field (for a simple introduction 
to this subject see reference [ll|). We see that the addition of matter with constants hi ~ 1 and /12 ~ 1 changes 
"slightly" the original geodesic equations for the Godel metric. We still have closed curved that are similar to the 
original Godel CTCs. Therefore we can say that the Godel CTCs are structural stable under the inclusion of the 
special matter represented by the energy-momentum tensor (|60[) . 

V. CONCLUSIONS 

In this work we verify that the closed timelike curves in Godel spacetime are stable under linear perturbations. We 
also show that the Godel spacetime has a stable structure under a special class of deformations. We found explicit 
CTCs in the deformed spacetime and proved that closed timelike geodesies do not exist. The energy-momentum 
tensor of the deformed spacetime was studied in some detail, specially we examined the conditions to have exotic and 
usual matter. We studied the stability of the new CTCs under linear perturbation and found that these curves are 
also stable. We tested these curves in a spacetime with exotic matter and find the same properties of stability as in 
the case of ordinary matter. We also find a kind of structural stability of the CTCs. 
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Godel in his seminal article [l[ mentions that his metric has five isometries [l|; Kundt 0] shows explicitly four out 
of the five Killing vectors of the above mentioned metric in Cartesian coordinates, 

C(o) = |, C(i) = §~ x , C(3) = §~ z , (81) 



^= X 0- X +y d-y- (82) 



The fifth Killing vector is not so trivial we find, 



c w = T (y " 1) ^ + 2 (1 - y + ^te +a *W (83) 

For a discussion of the topology of Godel metric and the five Killing vectors in a different system of coordinates, 
see [1^. 

The five Godel Killing vectors plus the discrete symmetry of reflection on the z — plane ( z — > — z) give us a 
family of metrics, that we will named the Godel family, 

ds 2 = kidt 2 + —dtdx + ^dx 2 + ^dy 2 + k 5 dz 3 ; (84) 

y y y 

where fej, i = 1, . . . , 4 satisfy the following relations: 

fc 2 - !L-ki = 0, ^jk 2 - k 3 + fc 4 = 0. (85) 

The corresponding constants that appear in the metric (|5^)) do not satisfy the above relation. In fact, only four 
of the above mentioned five Killings vectors (C(o)j •••■> C(3)) are symmetries of (1591) . therefore this last metric can be 
considered as a deformation of the Godel metric. 
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